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We propose technique of simultaneous excitation of disordered qubits providing an effective sup¬ 
pression of inhomogeneous broadening in their spectral density. The technique is based on applying 
of optimally chosen 7r-pulse with non-rectangular smooth shape. We study excitation dynamics of 
an off-resonant qubit subjected to strong classical electromagnetic driving field with the fast ref¬ 
erence frequency and slow envelope. Within this solution we optimize the envelope to achieve a 
preassigned accuracy in qubit synchronization. 


I. INTRODUCTION 

The investigation of qubit ensembles reveals analogies 
with quantum optics effects [1-3] and possibilities for 
construction of quantum computers and simulators [4— 
]. Solid state realizations of qubit ensembles are super¬ 
conducting Josephson circuits [4, 8 , 9], nitrogen-vacancy 
(NV) centers in diamond samples [10-12], or nuclear 
and electron spins realized as 31 P donors in 28 Si crys¬ 
tals [13] and Cr 3+ spins in AI 2 O 3 [14]. The coupling of 
qubit ensembles with a superconducting microwave res¬ 
onators results in the formation of sub-wavelength quan¬ 
tum metamaterials [15—19]. The long-range interaction 
through a photon mode results in the formation of collec¬ 
tive qubit states in such nretanraterials [20, 21], as Dicke 
model describes. One of crucial distinctions of artifi¬ 
cial qubits from natural atoms is that their excitation 
energies are in many cases tunable in situ by external 
magnetic fields. Beside of the tunability, another prop¬ 
erty is a disorder in excitation frequencies and, as a con¬ 
sequence, inhomogeneous broadening of the density of 
states in qubit ensembles. This is related to fundamental 
mechanisms such as exponential dependence of excita¬ 
tion energy on Josephson and charging energies in super¬ 
conducting qubits or spatial fluctuations of background 
magnetic moments [1 ] in systems with NV-centers. 

Disordered spectrum of collective modes offers multi- 
mode quantum memory, where information about pho¬ 
ton state is encoded as a tunable collective qubit mode 
[23, 24]. The storage and retrieval protocols were pro¬ 
posed in Refs. [25-27] and based on spin-refocusing tech¬ 
niques [28, 29] or successive magnetic field gradients. In 
the context of quantum memory the unavoidable spectral 
broadening in qubit excitation frequencies provides mul¬ 
timode performance from one side, but from the other 
side this is one of limiting factors affecting coherence 
times. Therefore, the development of techniques of ef¬ 
fective suppression of the disorder in qubit frequencies 
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and synchronization of their dynamics is an important 
problem. For instance, one of the options is atomic 
frequency comb (AFC) technique which could be ap¬ 
plied to rare-earth-metal-ion qubit ensembles [30]. This 
method is based on frequency-selective optical pumping 
and subsequent transitions to metastable auxiliary hy- 
perfine states. Another way of solution of this problem 
was demonstrated in Ref. [10] as ’cavity protection’ ef¬ 
fect in NV-centers. The effect is related to a decreasing 
of a relaxation rate of collective qubit modes which is 
proportional to the spectral broadening. 

Our research is inspired by one of the key ideas of Ref. 
[ 10 ]: the succession of microwave rectangular pulses can 
serve as efficient method for excitation of disordered NV- 
centers from the ground to the excited state. In our pa¬ 
per we study the possibility of simultaneous qubit exci¬ 
tation by a single non-rectangular 7r-pulse, rather than 
the sequence mentioned above. We observe that the op¬ 
timized non-rectangular shape of 7r-pulse provides an ef¬ 
ficient tool for suppression of the disorder effects as well. 
It allows to excite qubits within a wide detuning range 
with almost 100% probability. In contrast to AFC meth¬ 
ods, this technique does not require using of auxiliary 
levels transitions. 

We assume that the 7 r-pulse is realized as electromag¬ 
netic signal f(t)e~' lult of a carrying frequency ui being al¬ 
most in resonance with the qubit excitation frequency. In 
our solution we perform the optimization of the envelope 
shape /(f) in an experimentally relevant class of smooth 
functions, which guarantees that higher energy levels of a 
qubit are not affected. We expect that this technique can 
be applied to disordered systems with strong qubit-cavity 
couplings like NV-centers or superconducting metamate¬ 
rials, as well as to the atomic clock devices [31] as a tool 
for the preparation of a particular atomic state. 

II. DEFINITIONS 

We address to the possibility of simultaneous exci¬ 
tation of disordered qubit ensemble coupled to pho¬ 
ton transmission line being the source of the driving. 
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Qubits are assumed to be non-interacting with each 
other and long lived in comparison to 7r-pulse dura¬ 
tion time r <C t v . The absence of qubit-qubit inter¬ 
actions means that we can study dynamics of a single 
off-resonant driven qubit. We fix carrying frequency 
ui and assume that the qubit energy e can be var¬ 
ied reflecting the spectral broadening. Neglecting the 
qubit decoherence we solve the Schrodinger equation only 
idt\ip) = (H q +H ext )\ip), where unperturbed qubit Hamil¬ 
tonian is H q = e(<To + <j z )/2 and the external driving is 
H ex t = (, f{t)e~ lut a + + f*(t)e lult a-)/ 2. We define wave 
function of the qubit state in w-rotating frame as 


i m) = 


V f3(t)e~ iSt ^ 2 


where the detuning frequency is S = e — ui. The Hamil¬ 
tonian H of the driven qubit in this rotating frame reads 





(1) 


We assume that the qubit is close to cavity resonance e ~ 
ui and consider the evolution of the qubit wave function 
within the time interval 0 < t < r starting from the 
ground state |^>( 0 )) = | g) at the initial moment of time 
t = 0. The 7 r-pulse time r is considered as a fixed value. 
At this point we define frequency 



T 


which is the main scale in our consideration along with 
the detuning S. The frequency H has the transparent 
physical meaning: this is frequency of Rabi oscillations of 
the resonant qubit with <5 = 0 under the constant driving 
amplitude given by F a (t) = Q,e~' lbJt . Hence, time r is 
the half of the Rabi period associated with rectangular 
7 r-pulse F o (0 < t < t) exciting the system from | g) to |e). 
Non-zero detuning, related to inhomogeneous broadening 
or spread in qubit frequencies, does not allow to achieve 
full qubit excitation if envelope shape f(t) is constant. 
In the following consideration we modify f(t) at the time 
interval 0 < t < t to more complicated non-rectangular 
shape f(t) ^ const to achieve higher efficiency in near- 
to-resonance qubit excitation. 

Schrodinger equation with the Hamiltonian (1) al¬ 
lows analytical solutions only in several particular cases. 
The basic one is the constant driving amplitude f(t) = 
f = const and arbitrary detuning S which corresponds 
to damped Rabi oscillations of the frequency flu = 
y/ f 2 + 6 2 . In this case the evolution of the wave function 
being in the ground state at t = 0 reads 

/ -^sinf2 fl t/2 \ 

I m) = (2) 

\cos H^t/2 + sin flat/2J 

One can see that detuning reduces the maximum of ex¬ 
citation probability. This effect results in impossibility 


of synchronization of qubit ensemble by the rectangular 
shape of the driving envelope. 

Our further consideration is based on another exact 
solution, which holds for on-resonance driving regime S = 
0 and arbitrary real valued f(t). The time evolution of 
the ground state within this solution reads 




( • . v(t)\ 

—i sm- 

2 

<p(t) 


V 


cos 


7 


( 3 ) 


where the phase tp(t) is given by the time integral 


t 

¥>{t) = J 

0 

The 7 r-pulse condition, which is the inversion of a qubit 
occupation number, for this resonant case holds for 


ip(r) = 7r. (4) 

The last constraint (4) provides the class of real valued 
functions f(t) we are addressed to in the optimization 
procedures below. 


III. PERTURBATIVE SOLUTION 


As it was mentioned above, the exact solution is not 
known for arbitrary f(t) and non-zero detuning 5 ^ 0. 
Hence, we develop a perturbation theory based on treat¬ 
ing the (fcr-yQ-terms in the Hamiltonian (1) as small per¬ 
turbation and considering the exact solution (3) at S = 0 
as the zero order approximation. We end up with the 
following recursive equations forming the perturbation 
theory by <5 

(5) 

( 6 ) 

The full solution reads 

W) cos _ i y (t) S in Ah 

W)/ " ^((IcosAl-ixWsin^J 

Assuming that the qubit was in the ground state at the 
initial moment of time |^>( 0 )) = | g), the solution for x(t) 
and y(t) is given by 5 n series with nested integrals as the 
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coefficients 


t 

S f 

x(t) = - - / dt isin^(ti) (8) 

o 

+ J dti J dt 2 sin [<p(t 2 ) - <p(*i)] + • • • 

o o 

t 

St f 

y(t)= 1 + — / dticos^(ti) (9) 

o 

2 t 1 1 

- -4 -J dtl J dt<2 cos + • • • 


0 0 


These are general equations of the perturbation theory 
we use to find optimal shape of a 7r-pulse /(f) provid¬ 
ing simultaneous qubit excitations. Assuming envelope 
of the driving as smooth function at t = 0 and t = r 
we model /(f) as superposition of finite number of N 
sine functions, where O comes as a factor both in sine 
arguments and the driving amplitude 


N 

m = n E k 2n _i sin(?z — l/2)Of. (10) 

n =1 


According to (7) the wave function after the 7r-pulse 
takes the following form 


/a(r)\ /— iy{r)\ 


( 11 ) 


We optimize numerically a finite set of coefficients k 2 n-i 
and require that ground state amplitude at t = r is zero 
up to S N in (5-expansion, i.e. /3(r) = 0(S N ). It is pos¬ 
sible if two conditions are fulfilled: (i) resonant qubit is 
excited, i.e. ip(r) = i r, expressed in terms of k 2n -i as 


^ K2n-1 

^ 2 n- 1 

n—1 


= 7T, 


( 12 ) 


and (ii) off-resonant qubit excitation almost does not de¬ 
pend on detuning, i.e. x(r) = 0(S N ) in (11). The re¬ 
quirement x{t) = 0(5 N ) can be reduced to the set of 
N—1 equations corresponding to vanishing of <5" terms 
at n < N — 1 in the expansion (8), if we represent it as 
x(t,S) = X) c n8 n , where 

n—1 


T 

ci = J dtismip(ti) = 0 

o 

c 2 = j dti J dt 2 sin [tp(t 2 ) - v{t\)\ = 0 


T 1 1 


0 0 


(13) 

(14) 


To summarize, our perturbative solution consists of 
finite system of N equations (12, 13, 14, ...) providing 
smooth solution for /(f). The precision order of the tech¬ 
nique is given by N and ensures that the excitation prob¬ 
ability 71 |(t, <5) of an off-resonant qubit is close to unity 
up to small correction, namely n|(r, 6) = 1 — |x(r, <5)| 2 . 
This correction is nothing but the precision of the tech¬ 
nique and is given by the probability of the ground state 
qubit occupation tio(t,S) = \x(t,S)\ 2 . 


IV. RESULTS 

A. 7r-pulse optimization scheme at N = 3 

In this subsection we provide numerical results for the 
optimization of 7 r-pulse constructed from N = 3 terms 

/(f) = II(fci sinffI/2 + sin 3tfl/2 + k$ sin 5ff2/2). 

We start from the numerical solution for k% and k§ using 
(13,14) and after that we choose k-\ according to 7 r-pulse 
condition constraint (12). In the fig.l we plot two sets of 
curves in coordinates (ks,k 5 ): (i) red curves correspond 
to vanishing of the linear in <5 term in x(t), i.e. Ci = 0, 
according to Eq. (13), hence, in this case the residual 
part for the ground state amplitude is x(t) = 0(<5 2 ); (ii) 
blue curves correspond to vanishing of the quadratic term 
c 2 = 0 in x(t), given by Eq. (14). Note, that under the 
last condition c 2 = 0 the linear in S term may survive 
(ci 7 ^ 0). Crossing points of these blue and red set of 
curves satisfy both of the conditions (i, ii). These points 
correspond to synchronization of qubits excitations with 
the precision no(r, 6) oc (5/tt) 6 . The third parameter ki 
is found from (12) where we set N = 3 and take /c 3 , k$ in 
accordance with blue and red curves crossing points, e.g. 
that one is marked as greed dot. This green point corre¬ 
sponds to envelope /(f) having smallest maximum value 
which provides the effective suppression of the disorder. 


B. Synchronization of the qubits excitation vs 
detuning 

The higher order schemes are build straightforward 
around the above solution at N = 3. In this section we 
collect all the results for N = 1,2, 3,4 order schemes in 
the driving envelope function (10). In the left column of 
the fig. 2 we plot the optimized shapes of 7r-pulses found 
within the above perturbative approach for a given trun¬ 
cation number N. In the right column we show plots 
illustrating time evolution of qubit excitation dynamics 
nf(t) within 7r-pulse duration time 0 < t < r. Solid 
curves in the right column in fig. 2 correspond to reso¬ 
nant driving (5 = 0, while dashed ones are related to the 
detuning S = Q and 20. The dynamics starts from the 
ground state at t = 0 and grows significantly at the half of 
the 7r-pulse duration time r/2. The increase of the n\{t) 
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FIG. 1: Curves in (^’ 3 ,^ 5 ) plane formed by optimal 
values of sine amplitudes. Red curve corresponds to 
quadratic dependence of ground state amplitude on 
detuning. Crossing points of red and blue curves 
correspond to cubic dependence of ground state 
amplitude on detuning. Such a point nearest to zero is 
marked by green circle and studied in details. The value 
of fc-| is selected so that 7T-pulse is formed. 


resemble the response to singular driving at t = r/2, be¬ 
cause the limit of /(<) at N —► 00 corresponds to ideal 
7r-pulse with /Ar_>.oo(i) = n5(t — r/2), which is obviously 
not achievable experimentally. We stress that we work 
in the regime of finite N and amplitudes and treat the 
efficiency of this approach by means of deviation of the 
resulted n-j- (r) from unity with respect to non-zero detun¬ 
ing S. The last two plots in fig. 2 illustrate good efficiency 
of the corresponding 7r-pulse shapes: at N = 3 and 4 
we observe that the dashed curves are very close to the 
solid ones at t = t. This means that the inhomogeneous 
broadening is effectively suppressed in a Rabi frequency 
range oc Q and this technique offer the synchronization 
of qubit ensemble starting even from N = 3. 

In the fig. 3 we plot numerical results for the dependen¬ 
cies of ground state amplitude absolute values y/no(r) 
after 7r-pulse as a function of detuning S associated with 
a spectral broadening. It can be seen that the increase 
of th e schem e order N results in flattening of the curves 
for yj 77 . 0 (t ) around point <5 = 0. This flattening is quan¬ 
titative demonstration of the inhomogeneous broadening 
suppression. 

The fig. 4 is plotted in double logarithmic scale illus¬ 
trates the precision tiq{t) of the 7r-pulse technique pro¬ 
posed. This figure allows one to estimate the residual 
value of ground state amplitude for the given scheme or¬ 
der. For instance, for N = 3 the optimized 7r-pulse allows 
us to achieve the probability of qubit ensemble excitation 
up to n-f « 1 — 10 -3 at detuning values up to 0. 

The driving amplitude is limited in real experiment. 
To illustrate the effect of this limitation we plot the resid- 






FIG. 2: The shapes of driving amplitudes (on the left) 
and time dependence of excited level occupation 
number n^(t) (on the right) plotted for schemes of a 
different orders. The driving envelope f(t) shown in 
fig. 2a corresponds to non-optimized 7r-pulse and is 
plotted as a reference. 

ual ground state amplitude x(r) for the pulse of the same 
shape as mentioned above, but at different r so that max¬ 
imum value of the corresponding envelope does not ex¬ 
ceed Q. These curves are plotted as dashed lines in fig. 4. 

From the curves for y^ng (t) shown in fig. 4 we ex¬ 
tract the coefficients of power-law dependencies of resid¬ 
ual ground state amplitude for a given N. As one can see, 
the dimensionless combination no(6)(5/Cl)~ N does not 
depend on S. Dependence on N of logarithmically scaled 
value of this combination can be fitted by a straight line, 
as shown in fig. 5 

/ c \ 2 N 

n 0 (<5) = l —34(0.16-J (15) 

Eq. (15) is one of central results which show quantitative 
dependence of the precision on the order N and detuning. 
The small scaling factor 0.16 for S/£l shows that this 
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FIG. 3: Dependence of ground state amplitude after 
optimized 7r-pulse as a function of detuning in schemes 
of different orders N. 



FIG. 5: Ground state occupation normalized by 
appropriate power of detuning after optimized 7r-pulse 
as function of scheme order N involved in f{t). 


V. DISCUSSION 


technique based on sine representation (10) could provide 
the synchronization of ensemble even if the driving signal 
amplitude is less than a broadening by one order of the 
magnitude. 



FIG. 4: Dependence of residual ground state occupation 
number after optimized 7r-pulse on detuning value 
plotted in double logarithmic scale. Solid lines 
correspond to optimized 7r-pulses of the same duration 
r = 7r/fl. Dashed lines correspond to normalized 
envelopes having maximum absolute value equals to 
unity. In the latter case the 7r-pulse time is not fixed 
and increases with the increase of scheme order. 


The efficient excitation of NV-centers in diamond re¬ 
ported in Ref. [10] were achieved through the sequences 
of rectangular pulses with periodical switching of the 
amplitude sign. This technique demonstrates possibil¬ 
ity of excitation of strongly off-resonant qubits by a 
weak driving signal. We studied an opposite regime 
when the single non-rectangular 7r-pulse effectively sup¬ 
presses the disorder. The mechanism we are addressed 
to is related to synchronous excitation dynamics of two 
level system under a particular non-rectangular envelope 
shape f(t) of the 7r-pulse. We considered non-rectangular 
smooth shape of the 7r-pulse given by external electro¬ 
magnetic driving F(0 < t < r) = f(t)e~ lult with the 
envelope representable as the sum of N sine functions 
/(*) = k2n+i sin(2n + l)7ri/(2r) where the am¬ 

plitude and pulse duration are locked with each other 
r = 7 r/fi. The off-resonant response of a qubit to a non- 
rectangular signal can not be calculated exactly and we 
found the perturbative solution. Within this solution we 
proposed the method based on optimization of the set 
of N parameters fc 2 n+i which provide synchronous ex¬ 
citation of the off-resonant qubits. Note, that this op¬ 
timization is not direct expansion in sine basis of ideal 
7r-pulse in delta-functional form. The precision of this 
method, expressed in terms of qubit excitation number 
n-j-, is controlled by the order of the scheme and propor¬ 
tional to (5/Q) 2N . This scheme is efficient for the qubit 
energies falling into tunable spectral range estimated as 
the driving amplitude strength D. Within our solution 
we demonstrated that the 7r-pulse formed by N = 4 sine 
functions shows simultaneous excitation of qubits with 
the probability up to ss 1 — 10“ 5 for qubit frequencies 
ranging in ss u ± Cl. 

The sine expansion we used in this approach was based 
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on the experimental requirement of continuity of /(f) at 
initial and finite moments of time. Optimal envelope 
shape /(f) can also be found in another basis, say co¬ 
sine series or rectangular-based blocks. Our calculations 
shows that in these cases the results will be qualitatively 
the same as described above. Thus, the proposed method 
is quite general and can be tuned to meet requirement 
and restrictions of a particular experiment. 

To conclude, we propose the model of smooth shaped 
single 7r-pulse which can be applied to realistic disordered 
qubit ensembles coupled to a transmission line. Such a 
7r-pulse provides an effective suppression of the inhomo¬ 
geneous broadening and can be used as qubit synchro¬ 


nization technique. We expect that our findings could 
serve as a complementary methods to those reported in 
Ref. [10], where the sequences of rectangular pulses were 
used to increase the efficiency in the excitation of qubits 
within certain frequency range. We also mention that 
similar technique can be effectively applied to create, say, 
7r/2-pulse to prepare entangled states of the inhomoge- 
neously broadened qubits. 
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